Digital holography with ultimate sensitivity 
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We propose a variant of the heterodyne holography scheme, which combines the properties of 
off-axis and phase-shifting holography. This scheme makes it possible filter-off numerically the 
zero order image alias, and the technical noise of the reference. It is then possible to record and 
reconstruct holographic images at an extremely low signal level. We show experimentally that the 
sensitivity of the method is limited only by the quantum nature of photons. 



Digital holography is a fast-growing research field that 
has drawn increasing attention The main advantage 
of digital holography is that, contrarily to holography 
with photographic plates [|J , the holograms are recorded 
by a CCD and the image is digitally reconstructed by a 
computer, avoiding photographic processing [3|. Off- axis 
holography [4| is the oldest and the simplest configuration 
adapted to digital holography [1, H, 0] • 

In off-axis digital holography, as well as in photo- 
graphic plates holography, the reference or local oscillator 
(LO) beam is angularly tilted with respect to the object 
observation axis. It is then possible to record, with a 
single hologram, the two quadratures of the object com- 
plex field. However, the object field of view is reduced, 
since one must avoid the overlapping of the image with 
the conjugate image alias P. 

Phase shifting digital holography [3,[l3] records several 
images with different phase for the LO beam. It is then 
possible to obtain the two quadratures of the field in an 
on axis configuration even though the conjugate image 
alias and the true image overlap, because aliases can be 
removed by taking images differences. 

In this paper, we propose a digital holography tech- 
nique that combines off-axis geometry (introduced by 
Schnars et al. @), with the use of a sequence of images 
obtained with different phase shifts of the LO beam to 
record the hologram in amplitude and phase (as proposed 
by Yamaguchi et al. @). To get precise phase shift, Le 
Clerc's et al. [13] heterodyne technique is used. Using 
a spatial filtering method (Cuche's et al. @), the zero 
order image, and the noise attached to it, is filtered-off 
numerically. As shown experimentally, this combination 
of techniques makes it possible to record and reconstruct 
holographic images at a very low level of signal: 1 photo 
electron of signal per reconstructed image pixel during a 
whole sequence of 12 images (~ 1 s). This corresponds 
to the ultimate quantum limit. 

The setup is shown on FigJTJ It is similar to the one 
used in [13, IH . The main laser L is a Sanyo DL-7140-201 
diode laser (A = 785 nm, 50 mW for 95 mA of current). 
It is split into an illumination beam (frequency ojl , com- 
plex field El), and in a LO beam (u>lOj Elo)- The 
object we want to image is an USAF target in trans- 
mission, which is back illuminated. The object signal 
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FIG. 1: USAF target digital holography setup. L: main laser; 
BS: Beam splitter; AOM1 and AOM2: acousto optic modu- 
lators; BE: beam expander; M: mirror; A: light attenuator. 
USAF: transmission USAF target. CCD : CCD camera. 



field E is not shifted in frequency {u>l)- A set of opti- 
cal attenuators A (grey neutral filter) is used to reduce 
the illumination. The CCD camera (PCO Pixelfly digital 
camera: 12 bit, frame rate ujccd — 12.5 Hz, acquisition 
time T = 1/12.5 = 80 ms, with 1280 x 1024 pixels of 
6.7 x 6.7 nm) records the hologram of the object, i.e. the 
object (E) and LO (Elo) field interference pattern. Us- 
ing AOM1 and AOM2 acousto optic modulators (Crystal 
Technology: ^01/1,2 — 80 MHz), the optical frequency 
ojlo °f the LO beam, can be freely adjusted. To make 4- 
phases detection of the object field E, the LO frequency 
is adjusted to be ojlo = + ^>ccd/^ [I3|- Moreover 
the LO beam is tilted (angle 9 ~ 1°) with respect to 
the camera to object observation axis. Our holographic 
setup thus works both in off-axis and in phase shifting 
mode. 

A sequence of 12 CCD images Iq to In (measurement 
time 0.96 s) is recorded. Since the LO beam is phase 
shifted by tt/2 between two consecutive images, the ob- 
ject complex hologram H is obtained by summing the 
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CCD images with the appropriate phase shift: 



H 



(1) 



where m is the image index and j 2 = — 1. 

We have reconstruct the images of the USAF target by 
using the standard convolution method [1, [13] that yields 
a calculation grid equal to the pixel size. To calculate the 
convolution product, we have used the Fourier method, 
like in (ic| . To make faster Fourier calculation the 1028 x 
1024 CCD data are truncated to a 1024 x 1024 2D matrix. 

The reconstruction algorithm is the following. The real 
space hologram H (x, y, z) in the z = CCD plane is 
calculated by Eq[TJ The hologram in the CCD reciprocal 
plane (i.e. in the z = k-space) is obtained by Fourier 
transformation: 



H(k x , k y , 0) = FFT [H(x, y, 0)] 



(2) 



The k-space hologram at any distance z from the CCD 
is then: 



H[k x , ky, z) — K(k x , ky , z)H(k x , ky, z) 



(3) 



where K is the k-space kernel function that describe the 
propagation from to z. 



K(k x ,k y ,z) = e»*(*-"+*» a )/* 



where k — 2it/\ is the optical wave vector. The recon- 
structed image, which is the hologram in the object plane 
(z = D), is then obtained by reverse Fourier transforma- 
tion: 



H(x,y,z)=FFT~ 1 \H(k x ,k y ,z 



(5) 



FigGk shows the intensity of the CCD plane complex 
hologram (i.e. \H(x, y, 0)| 2 ). The USAF target is seen, 
but is blurred. FigOs shows the intensity of the CCD 
plane k-space hologram (i.e. \H{k x , k y , 0)| 2 ). The bright 
zone in the left hand side of the image corresponds to 
the true holographic image i.e. to the interference of 
the object field with the LO field. The zero-order and 
twin image Q alias contributions are very low. The zero- 
order alias, in the center of the k-space hologram (pixel 
512, 512), is weak, because the contribution of the LO in- 
tensity cancels out by making difference of images (EqfT]). 
Since the phase shift provided by the acousto optic mod- 
ulators is very close to a multiple of tt/2, the twin-image 
is also very weak. It is barely not visible on FigOj, al- 
though the display is in grey logarithmic scale. 

To select the relevant first order image information, 
and to fully suppress the zero-order and twin image 
aliases, we have used the k-space filtering method de- 
veloped by Cuche et al. We have selected, in the 
k-space 1024 x 1024 matrix, a 200 x 200 region of interest 
centered on the true image bright zone. Note that this se- 
lection is made possible by the off-axis geometry that has 




FIG. 2: Reconstruction of the USAF target image, (a) In- 
tensity image of the CCD complex hologram \H(x, y, 0)| 2 in 
linear grey scale, (b) Intensity image of the k-space hologram 
\H(k x , k y , 0)| 2 in logarithmic grey scale, (c) Intensity image 
of the truncated k-space hologram in logarithmic grey scale, 
(d) USAF target reconstructed image, i.e. intensity image of 



(4) the real space hologram in the object plane \H(x, y, D)\ 2 



linear grey scale. Images are 1024 x 1024 (a,b), or in 256 x 256 
(c,d) pixels. 



translated the true image in the left hand side of k-space 
matrix. The selected area is then copied in the center of 
a 256 x 256 zero matrix (zero padding). The calculation 
of the z = D k-space and real space holograms (Eq[3] 
and Eq[5]) are thus done on this 256 x 256 calculation 
grid. Figl^t shows the object plane k-space hologram (i.e. 
H(k x ,k y , z)) in the 256 x 256 grid. Fig^U shows the ob- 
ject plane real space hologram (i.e. H{x,y,z)) obtained 
by the Eq[5] reverse FFT, which is the reconstructed im- 
age of the USAF target. Note that the translation of the 
selected zone in the center of the k-space in FigJ5J; moves 
the reconstructed USAF target to the center of the image 
as seen in Fig(5Ji. 

We have assessed the sensitivity limit of this off-axis 
heterodyne holography technique by recording images of 
the USAF target at different levels of illumination. To 
get quantitative results, we have determined the absolute 
number of photo electrons that correspond to the signal 
beam. First, we have calibrated our set of attenuators 
(A) at the working wavelength (785 nm). These atten- 
uators were then used to change the illumination of the 
USAF target while the laser power remained constant 
throughout the experiment. For a high level of illumina- 
tion of the USAF target, and without LO beam, we have 
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FIG. 3: Reconstructed image of a USAF target in transmis- 
sion with low light illumination. Images are obtained with a 
signal field that corresponds to ~ 2.8 x 10 r (a), 2.8 x I0 5 (b), 
4 x I0 4 (c) and 6000 (d) photo electrons for the whole 12 CCD 
images (256 x 256 pixels displayed in linear grey scale). 



filter-off technical noise and to reach the quantum limit 
of photodetection. Note that it is still possible to per- 
ceive the USAF target on FigOi with only 0.15 photo 
electron per pixel (SNR~ 0.15). 



This result can be understood quite simply. In the 
case of a very weak signal, the majority of the technical 
noise is within the LO beam (because the signal beam 
is much weaker than the LO beam). This noise can be 
eliminated by a double filtering process. First, most of 
the noise is removed by taking the difference of succes- 
sive images (EqQ}. This is the phase-shifting filtering 
that occurs in the temporal frequency domain. But this 
method doesn't remove all the technical noise. Since the 
LO beam extends within a few spatial modes (i.e. a few 
pixels of the k-space), its noise is located near the center 
of the k-space. This means that the LO technical noise 
is flat field (it varies slowly in space). A spatial filtering 
similar to the one introduced by Cuche et al. [8j removes 
the rest of the technical noise which lies in the center of 
the k-space. This is the off-axis filtering that occurs in 
spatial frequency domain. 



measured the signal beam in photo electrons units: we 
have multiplied the CCD signal in digital counts (0...4095 
DC) by the CCD gain (G = 2.4 photo electrons per DC: 
PCO calibration). For other levels of illumination (for 
a given set of attenuators in the object arm), we have 
calculated the signal beam level in photo electrons units, 
from the attenuators calibration. 

Fig[3]shows the reconstructed images that are obtained 
at various levels of illumination. Figs|3] a to d are ob- 
tained with ~ 2.8 x 10 7 (a), 2.8 x 10 5 (b), 4 x 10 4 (c) 
and 6000 (d) photo electrons of signal for the whole set 
of 12 CCD images. It should be noticed that the total 
number of photo electron is not the relevant parameter to 
characterize the image signal. Since the reconstruction 
is performed with a truncation over a 200 x 200 pixels 
region of the k-space, the reconstructed image has about 
4 x 10 4 resolved pixels. The signal can thus be character- 
ized by the number of photo electrons per resolved pixel. 
Therefore, the images of Fig|3] correspond to 700 (a), 7 
(b), 1 (c) and 0.15 (d) photo electrons per pixel. For 
comparison, the image of Figj2ji corresponds to 3.5 x 10 4 
photo electrons per pixel. 

One photo electron per resolved pixel is the quantum 
limit of visibility of an image (i.e. the limit that is related 
to the quantum character of photons, and photo elec- 
trons). Figl3j: image corresponds exactly to this limit, 
where the SNR of the image is about 1. This demon- 
strate that our off-axis heterodyne technique is able to 



The digital holography technique presented here relies 
on the combination of off-axis and phase-shifting config- 
urations to record holograms with ultimate sensitivity. 
This technique lets one filter-off the local oscillator tech- 
nical noise and reach the quantum sensitivity limit. As 
demonstrated, we were able to make an image with a sig- 
nal of one photoelectron per pixel during a sequence of 
12 recorded camera frames. Although demonstrated in 
through transmission this technique is expected to work 
as well in reflection with diffusely scattering surface. It 
might be used to perform holography with at extremely 
low signal level, e.g. nano-object imaging... We have 
used a variant of this technique to detect the so called 
"tagged photons" in acousto-optic imaging [l3|], or the 
photons that are transmitted through the breast, in vivo 
|14| . With a high quantum efficiency camera, the tech- 
nique could also be used to perform quantum optics tests 
(non classical photon statistics, squeezed states). 
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